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Executive Summary 



Tlie Nortliem Region of the U.S. Forest Service 
(USPS Region 1) provides habitat for a number 
of land mollusk species that are endemic to tlie 
Nortliem Rocky Mountains, as well as several 
otliers diat aie restricted to tiie Pacific Nortiiwest 
wiUi disjunct populations in nortliem Idaho and 
nordiwestem Montana. Tlie goal of this inventory 
project, which began in 2005, was to fill species 
distribution gaps, test survey methods, and collect 
data for development of predictive distiibution 
models. 

Using tlie NatureSeive website and prioi" reports 
as a stalling point, we compiled a hst of 29 land 
snail taxawitliin tlie U.S. Porest Service Nortliem 
Region that weie globally ranked in 2005 as Gl- 
G3 or T1-T3, thereby meeting USPS Species of 
Concern (SOC) criteria. Two additional snail taxa 
weie included with global ranks of G5 and state 
ranks S1-S2, tliereby meeting criteria for USPS 
Species of Inteiest (SOI). We also compiled a list 
of eight slug taxa ranked in 2005 as G1-G3, and 
three additional slug taxa ranked G4-G5 andSl-S2, 
again meeting the respective USPS criteria for SOC 
or SOI. Nmeteen of these SOC/SOI species (11 
snails, 8 slugs) have been documented in Montana 
as of 2007. 

Montana Natural Heritage Program surveys in 
2005 included lands in botii Idaho and Montana; 
in 2006 and 2007 tliey were restiicted to Montana. 
Tliis report presents results of tlie 2007 surveys 
which focused on gathering additional distiibution 
data to aid in the development of predictive 
distribution models and maps. We conducted a 
total of 34 site surveys on National Porest units in 
Montana, primaiily targeting aieas with highest 
expected diversity of SOC and SOI taxa but also 
areas wheie prior surveys were lacking; these 
were mostly (but not exclusively) west of tlie 
Continental Divide. Site surveys weie distributed 
on the Montana Forests as follows: Kootenai (16), 
Lewis & Claik (8), and Lolo (10). 

We found SOC and/or SOI taxa at 55.9% (19) 
of the sites surveyed, and documented a total of 
38 locations for nine SOC taxa and two SOI taxa 



during our 2007 surveys. Tliese taxa included: 
Alpine Mountainsnail Oreohelix aJpina (1 site), 
Carinate Mountainsnail Oreohelix elrodi (1 site). 
Humped Coin Folygyrelio polygyrella (1 site), 
Pir Pinwheel Radiodiscus abietum (6 sites). 
Pale Jumping-slug HemphiUia camelus (4 sites), 
Maibled Jumping-slug HemphiUia danielsi (1 
site). Magnum Mantleslug Magmpeha mycopitaga 
(4 sites). Pygmy Slug Kootenaia burkei (9 sites). 
Smoky Taildioppei" Frvphysaott humile (6 sites), 
Lyi'e Mantleslug Udosarx Jyrata (1 site), and 
Sheathed Slug Zacoleus idahoensis (4 sites) . Most 
were found west of the Continental Divide in niesic 
forest habitats (e.g., western redcedar, western 
hemlock, mesic Douglas-fir, grand fir), altliough 
tlie two species of mountainsnail (Oreohelix) were 
fiom exposed rocky teiiain above treeline on tiie 
Scapegoat Plateau just east of tlie Continental 
Divide. We also collected distribution data on 19 
additional non-SOC/SOI species as we encountered 
tliem duiing oui" surveys. 

We collected new locations for two slug SOCs tliat 
were found in Montana for the first time during our 
2005 surveys: Pale Jumping-slug and Pygmy Slug. 
Our 2007 surveys also added several new Montana 
locations for a tiiird slug SOC, Smoky Taildiopper, 
which was documented in Montana only once prior 
to 2004. Several SOC snails and slugs were found 
at new locations during tiie 2005-2006 suiTeys, 
and five Global and State Ranks were reduced as 
a result; additional rank changes may occur as a 
result of tiie 2007 surveys. 

Predicted species distribution models were 
developed at 90 m pixel resolution, using our point 
location data and 11 environmental data layers 
chosen for theii' assumed relationship to land snail 
requirements. Model output indicates tiiat the 
greatest number of gaps in known species ranges 
(especially Species of Concern) may be filled in 
witii additional surveys focused on portions of 
nortiiwestein Montana, the BitteiToot Mountains 
in Mineral and Ravalli Counties, Beaverhead, 
Madison, and Silver Bow Counties, and portions 
of the island mountain ranges in eastern Montana. 
In addition to tiiese gaps, a summary of elevations 



ui 



on U.S. Forest Seivice lands in Montana and 
elevations of mollusk survey locations in 2005- 
2007 indicates tliat high elevation sites have 
been under-sampled. We intend to field-test tlie 
predictive distribution models with surveys focused 
on these regions during the 2008 field season. 
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Introduction 



Tlie Nortliem Region encompasses tlie ranges 
of a large number of land mollusk species 
endemic to the Nortliem Rocky Mountains, and 
several additional species restricted to tlie Pacific 
Nortliwest witli disjunct populations in nortliem 
Idaho and noitliwestem Montana (Frest and 
Johannes 1995, 1997, 2001; Hendiicks 2003; 
Hendi'icks et al. 2006). Tlie US Forest Service is 
required under the National Forest Management 
Act (1976) and Code of Federal Regulations (CFR 
1985) to maintain a diversity of plant and animal 
species. Inventory is a fiist step in the evaluation 
of landscapes and tlieir likelihood of supporting 
populations of animal species of conservation 
concern. Pursuant with tliis legislation and 



associated regulations, the Nortliem Region 
initiated surveys in 2005 for a suite of land 
mollusks listed as Species of Concern in Montana 
andldaho (Hendricks et al. 2006). 

Objectives of the 2005 inventory included 
filling species distribution gaps, testing survey 
methodology, and collecting geospatial and 
habitat data for the development of predictive 
habitat and distribution models that can aid fiiture 
survey efforts. The survey was continued in 2006 
and 2007 with tlie same objectives, and witli the 
Montana Natural Heritage Program restricting 
its field effort to the nine National Forests within 
Montana {Hendiicks et al. 2007). 



Methods 



Prior to conducting field surveys in 2005, we 
searched tlie pnblislied aiid gray literature to 
compile a list of high-priority ''target" species 
(globally aiid state rare species in Montana, and 
globally raie species in Idaho). Primaiy sources we 
used for this compilation included Pilsbry (1939, 
1948), Fiest and Johannes (1995, 1997, 2001). 
and Hendi'icks (2003). Tliis resulted in a list of 
41 species and subspecies (31 snails, 10 slugs) we 
considered to be of conservation concern (Table 
1 ill Hendi'icks et al. 2006); 12 of these taxa are 
ranked less than G3 (see Appendix A for ranks 
definitions). We tlien generated a list of general 
habitat associations for die high-priority species 
(Table 2 in Hendricks et al. 2006), to help us 
prioiitize habitats for our suiveys during 2005 - 
2007. Exti'emely limited infomiation for mollusks 
east of the Continental Divide in Montana made 
this process more problematic for tlie high-priority 
species tliat occur or might occur in that region. 
Of the 41 species and subspecies identified as of 
conservation interest on or near tlie USPS Nortliern 
Region, 19 (11 snails, 8 slugs) occur in Montana 
and were the primary focus of survey efforts in 
2006-2007. 

In 2007, we conducted field suiTeys for land 
mollusks in late July during a high-elevation 
mountain tiaverse, and in September to eai'ly 
October, when tlie weatlier was most suitable at 
lower elevations (cool and moist) for finding active 
snails and slugs. In 2007, unlike prior yeais, we 
concentrated our suiveys on filling geographic 
data gaps, especially tai'geting ai'eas along or 
neai' political boundaiies. We also restiicted our 
attention to ai'eas near or west of the Continental 
Divide, because tliis is the region with the largest 
number of Species of Concern (SOC) . Tlius survey 
effoit in 2007 was non-random and concentiated 
in tluee ai'eas within tlie Noitliem Region: 1) the 
nortliern portion of tlie Kootenai National Forest 
neai" tlie Britisli Columbia border, 2) portions of 



tlie Lolo National Forest neai' the Idaho bordei\ and 
3) a backcountry section of the Lewis and Clark 
National Forest in the Sc^egoat Wilderness. 

We selected sites for surveys based primaiily on tlie 
presence of perennial water, moist mature conifer 
forest, aspen, and/or limestone talus or otlier rock 
outcrops. At each site, we conducted timed surveys 
while searching under leaf litter, dead wood and 
bai'k, rocks imbedded in tlie ground, or digging into 
talus. Usually witliin a survey site we seaiched 
several locations with habitat features (such as 
biyophyte mats, dead wood and imbedded rocks, 
or talus slopes) considered by experts to be favored 
by snails and slugs, often concentrating seaiches in 
riparian zones. 

We recorded a vai'iety of habitat and site 
information at each survey location on standai'dized 
data forms (Appendix B), and digital photographs 
were taken at all sites. Survey data from 2007 
have been enteied into tlie Montana Natural 
Heritage Program Point Observation Database 
(POD); copies of the Idalio POD data collected 
in 2005 were sent to tlie Idalio Consei'vation Data 
Center (CDC) in Boise. We collected voucher 
specimens of all SOC taxa we discovered, as well 
as representatives of many other non-SOC taxa; 
vouchers were preserved in 95% ETOH in order to 
pemiit future genetic analyses. We sent SOC slug 
vouchers collected in 2005 to taxonomic experts. 
Tlieir identifications were verified, and we used tliis 
knowledge to make species determinations for die 
2006-2007 material. 

Point observation data was used to generate 
predictive distiibution models and map outputs. 
Tlie modeling effort and inteipretation of an 
example output aie included in Appendix E and 
outputs for individual species ai'e available on tlie 
Montana Field Guide website: http://fieldguide. 
nit.aov 



Results and Discussion 



Overview 

We conducted a total of 34 site surveys in 2007 
(Appendix C). Tliese were distributed on tlie 
Montana Forests as follows: Kootenai (16), Lewis 
& Clark (8), andLolo (10). Land mollusks were 
found at all sites surveyed, represented by one to 
ten species per site (30 species total). Of these, 
only thiee species (all non-SOC/SOI) were found 
on all tliree Forests: Forest Disc (Discus vhitneyi). 
Brown Hive {Encomdusfnhus), and Spruce Snail 
(MicrophysuJa ingersoUh). We documented a total 
of 23 species on the Kootenai National Forest, 
seven species on tlie Lewis and Claik National 
Forest, and 18 species on tlie Lolo National Forest. 
SOC and/or SOI taxa weie documented at 1 9 
(55.9%) of the sites, mostly west of tlie Continental 
Divide. Tlie aieas surveyed addressed some of tlie 
gaps identified in our report of the 2006 inventory 
(Hendricks etal. 2007). 

We documented 38 locations for nine USFS 
Region 1 SOC taxa and two SOI taxa during our 
2007 surveys (Table 1): Alpine Mountainsnail 
Oreohelix alpina (1 site), Caiinate Mountainsnail 
Oreohelix ehodi (1 site), Humped Coin Folygyrella 
polygyrelia (1 site). Fir Pinwheel Radiodiscits 
abietnm (6 sites), Pale Jumping-slug i/^w//?M//a 



camelus (4 sites). Marbled Jumping-slug 
HemphiUia daniehi (1 site), Magnum Mantleslug 
Mognipelta mycophaga (4 sites). Pygmy Slug 
Kootena/a bnrkei (9 sites). Smoky Taildioppei' 
Frophysaon humile (6 sites). Lyre Mantleslug 
Udosarx iyrata (1 site), and Sheatlied Slug 
Zacolens idahoensis (4 sites). Most locations are 
fiom west of tlie Continental Divide in niesic forest 
habitats (e.g., western redcedai", western hemlock, 
mesic Douglas-fir, grand fir), but locations for tlie 
two mountainsnails (Oreohelix) were exposed 
rocky limestone teirain above tieeline on the 
Scapegoat Plateau just east of tlie Continental 
Divide (Appendix D). 

In 2007, we collected additional location data for 
two SOC slug species discovered on tlie Kootenai 
National Forest in 2005 and new at that time to 
tlie Montana mollusk fauna. Pygmy Slug has 
now been documented at 20 sites; one of tlie new 
Pygmy Slug sites is on tlie Lolo National Forest, 
making four total sites on tlie Lolo and 16 on 
tlie Kootenai. Pale Jumping-slug has now been 
documented at seven sites, all on tlie Kootenai 
National Forest. We documented an additional 
location for the Marbled Jumping-slug on the Lolo 
National Forest in Mineral County. Tlie St. Regis 



Table I. Manber ofsiii-vey sites where Species of Concern 
Forests in Montana during the 2007 siu-vey (n = 34 sites). 


land mollusks were detected on Noithem 


Region 


SPECIES 


G 

RANK 


TOTAL 
SITES 


MONTA 


NA' 


KO 


LC 


LO 


SNAILS 


Alp ill e Mountain snail 


Oisohehx alpina 


Gl 


1 




1 




Carinate Mountainsnail 


Oi-eohehx ebvdi 


Gl 


1 




1 




Humped Coiu 


Polygyj-ella polygymlla 


G3 


1 






1 


Fii' Pinwlieel 


RadiodiscH^ abietiim 


G4 


6 


1 




5 


SLUGS 


Pale Jumping-slug 


HemphiUia camelus 


G4 


4 


4 






Maibled Jumping-slug 


He}nphiUia damehi 


G2G3 


1 






1 


Pygmy Slug 


Kootenaia bwkei 


Gl 


9 


8 




1 


Magnum Mautleslug 


Magmpelta mycophaga 


G3 


4 


3 




1 


Sniokey Taildi'opper 


Frophysaon humile 


G3 


6 


3 




3 


Lyre Mantleslug 


Udosarx iyrata 


G2 


1 






1 


Sheathed Slug 


Zacolens idahoensis 


G3G4 


4 


4 







' Montana Forcsls codes: Kooteuai (KO), Lewis & Clark (LC), Lolo (LO). 
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Figure I. Map oj noithwestein Montana show mg locations oj land moUiLsks during 2Q07. 



River appears to be the barrier separating this from 
tlie congeneric Pale Jumpiag-slug, but more field 
(and lab) workwould help detemiine tlie ranges of 
tliese two species in Montana. 

Tlie 2007 suiTey added six Montana locations 
to the 31 from 2006 for the Smoky Taildi'opper, 
tliereby bringing the total locations to about 41 
for a G3 slug which was documented in Montana 
only once prior to 2004. We also documented four 
new locations of Magnum Mantleslug. bringing 
tlie total Montana locations to about 25 (Hendi'icks 
et al. 2006). To date, the slug witli tlie fewest 
repoited localities in Montana (otiier tlian the 
recently-documented Pale Jumping-slug) is the 
Lyi'e Mantleslug, known from just seven sites, even 
tliough it was first documented in the state in 1965 
(Russell and Webb 1980); four of the sites were 
discovered during 2006-2007. 

As a result of the 2005 sui"veys, tiie Global Rank 
of Humped Coin changed from G2G3 to G3 {see 
Appendix Afor definitions), the Fir Pinwheel 
changed from G3 to G4, Pale Jumping-slug 
changed from G3G4 to G4, the Pygmy Slug 
changed from Gl G2 to G2, and Smoky Taildroppei' 
changed from G2 to G3. Updated Global and State 
Ranks are included in Appendix E. Table 1 . Tliese 
updates are a result of tlie 2005-2007 surveys, tlie 
detection probability and site occupancy sui^eys 
conducted in 2006 and 2007, and the predicted 
distribution models produced in 2007. In summaiy 
tlie 2005-2007 Northern Region surveys have made 
a significant contribution to our understanding of 
tlie cuirent status of several land mollusk species of 
conseiTation interest in Montana. 

Finally, we recorded 19 additional land mollusk 
species (including exotics) as we encountered 
tliem during our 2007 suiveys. These species ai'e 
not currently recognized as SOC or SOI, nor aie 
they likely to merit such status, and will not be 
discussed fiirther in this report. Distribution maps 
showing locations where we found tliese taxa 
can be viewed at the Montana Natural Heritage 
Program Tracker website http://nlip. mis. state. 
iiit.us/Trackei\ Available for viewing aie all of 
our suivey records for all species, including two 
of the Chrome Ambersnail {Catiftella i-ehderi) 



from Cai'bon and Fergus counties, Montana. 
Species of Catmelia are impossible to identify to 
species based on shells alone (T. Pearce personal 
communication), so our identification of shells 
from tliese sites remains tentative, and influenced 
by one prior Montana record from Meagher County 
(Pilsbiyl948). 

Species Accounts 

Alpine Mountainsnail (Oreohel/x alpuia): We 
found tliis species scattered among several 
locations at one site {tlie Scapegoat Plateau NW 
of Scapegoat Mountain) between 7790-8075 ft 
elevation in Lewis and Claik County, on the Rocky 
Mountain RD of tlie Lewis & Claik National Forest 
(Appendix C). Tliis is the easternmost documented 
location for tlie species so fai", an extension of about 
38 miles {61 km) ESE of tlie Rumble Lake site in 
tlie Swan Range (Fairbanks 2002). The Alpine 
Mountainsnail is a Montana endemic, known only 
from the Mission and Swan ranges prior to 2007 
(Elrod 1903b; Hendricks 1998; Fairbanks 2002), 
and tied to alpine and high-elevation subalpine 
limestone talus. Approximately 60 live individuals 
and 100 shells were found. No locations of 
sympatiy with tlie Carinate Mountainsnail were 
found on the Scapegoat Plateau. Additional details 
regai ding this species ai'e presented in Appendix D. 

Carinate Mountainsnail (Oreobelix eirodi): We 
found this species scattered among several 
locations at one site {tlie Scapegoat Plateau NW 
of Scapegoat Mountain) between 7850-7960 ft 
elevation in Lewis and Claik County, on the Rocky 
Mountain RD of tlie Lewis & Clark National Forest 
(Appendix C). Tliis is tlie easternmost documented 
location for the species so far, an extension of 
about 44 miles (71 km) ESE of tlie Dog Creek site 
discovered in the Swan Range in 2005 by Henning 
Stabins. Tlie Carinate Mountainsnail is a Montana 
endemic, known only from the Mission and Swan 
ranges prior to 2007 (Elrod 1903a; Fairbanks 1984; 
Hendj'icks 1998, 2003). In these ranges it is found 
primarily in ai'gillite talus witli or witliout overhead 
canopy, and well-below ti'eeline. On the Scapegoat 
Plateau it was found under alpine limestone 
talus and blocks witli no overhead canopy. 
Approximately 100 live individuals and 400 shells 
were found. No locations of sympatiy witli tlie 



Alpine Mountainsnail weie found on the Sc^egoat 
Plateau. Additional details regaiding this species 
are presented in Appendix D. 

Humped Coin (Polvgyrella polvgyreUaY We 
found this species at one site at 4100 ft elevation 
in Mineral County, on the Superior RD of tlie 
Lolo National Forest (Appendix C). The Humped 
Coin, first described from Montana and Idaho by 
Bland and Cooper (1861) and Cooper (1868), is 
also present in adjacent Washington and Oregon 
(Frest and Johannes 1995, 2001). All known 
Montana sites are clustered in Sanders and Mineral 
counties (Hendricks 2003, 2005: Hendiicks et al. 
2006, 2007). Canopy at tlie 2007 site included 
western redcedai", western hemlock, grand fir, and 
Engelmann spiiice. One live animal was found 
under a > 20 cm rock at this site. 

Fir Pinwheel (Radiodiscus abietum): We found 
this species at six sites between 3140-4960 ft 
elevation, in Lincoln, Mineral, Missoula, and 
Sanders counties, on UieLibbyRD of tlie Kootenai 
National Forest, and the Missoula, Ninemile, and 
Plains RDs of the Lolo National Forest (Appendix 
C). The Fir Pinwheel is restricted to northern 
Idaho, western Montana, and adjacent parts of 
Oregon and Washington (Biiinson and Russell 
1967; Frest and Johannes 1995, 2001; Hendricks 
2003, 2005; Hendricks et al. 2006, 2007). Canopy 
at the 2007 sites included western redcedar, grand 
fir, Douglas-fir, alder, black cottonwood, and 
western larch. Single live individuals were present 
at each site, mostly under downed wood, but also 
rocks and leaf litter. 

Pale Jumpins-slus jHemphiUia canielusy We 
found this species at four sites between 2980-4800 
ft elevation, in Lincoln County, on tlie Rexford and 
Tliree Rivers RDs of tlie Kootenai National Forest 
(Appendix C). Tliis species was first documented 
in Montana during tlie 2005 survey (Frest and 
Johannes 1995; Hendiicks 2003; Hendricks et al. 
2006, 2007). The Pale Jumping-slug appears to be 
restiicted to nortliem Idaho, and adjacent parts of 
Washington, British Columbia, Alberta, and now 
Montana (Frest and Johannes 1995, 2001; Forsyth 
2004). Frest and Johannes (1997, 2001) suggested 



tliat individuals from tlie Lower Salmon River 
drainage in Idaho might represent a taxon distinct 
from tliat to the nortli, but this possibility has not 
been resolved. Canopy at the 2007 sites included 
western redcedai", subalpine fir, Engelmann spruce, 
western hemlock, western larch, lodgepole pine, 
ponderosa pine, and aldei\ Between three and eight 
individuals (20 total) were found at each site under 
downed wood, leaf littei; and rock. 

Maibled Jumping-slug (HemphiUia danielsiY We 
found this species at one site, at 4100 ft elevation 
in Mineral County, on the Superior RD of tlie Lolo 
National Forest (Appendix C). This species was 
first documented in Montana in 1912 (Vanatta 
1914; Frest and Johannes 1995; Hendricks 2003). 
Until recendy, tlie global range was exchisively 
tlie Bitterroot Mountains. Tlie Maibled Jumping- 
slug appears to be restiicted to extieme western 
Montana soutli of the St. Regis River, neai' the 
state line witli Idaho (Frest and Johannes 1995, 
2001: Hendiicks et a!. 2006, 2007), altliough there 
is an older HemphiUia record near Flathead Lake 
tliat has been tentatively assigned to this species; 
tlie species also occurs in adjacent Idaho, but the 
distribution as yet is not well-defined. Canopy at 
tlie 2007 site included western redcedai; grand fir, 
Engelmann spiiice, and western hemlock. A single 
individual was found under downed rotten wood. 

Pygmy Slug (Kootertaia burkei) : We found 
tliis species at nine sites between 2680-4800 ft 
elevation in Lincoln and Mineral counties, on 
tlie Libby, Rexford, and Three Rivers RDs of 
the Kootenai National Forest and the Superior 
RD of tlie Lolo National Forest (Appendix C). 
Only recently was this species first discovered 
and described, fi^om five sites in nortliem Idaho 
(Leonard etal. 2003). It was documented in 
Montana for tlie first time at foui" sites during the 
2005 survey (Hendiicks et al. 2006), and at seven 
more sites in 2006 (Hendiicks et al. 2007). Canopy 
at the 2007 sites included western redcedai", grand 
fir, Engelmann spruce, western hemlock, subalpine 
fir, lodgepole pine, pondeiosa pine, alder, and 
western larch. Up to 26 individuals (41 total) were 
found on and under downed wood, under rock 
fiagments, and among leaf litter. 



Magnum Mantlcslug {Magnipelta ntycophagdy. 
We found tliis species at four sites between 3140- 
4500 ft elevation in Lincoln and Mineral counties, 
on the Fortine and Libby RDs of the Kootenai 
National Forest and tlie Superior RD of the Lolo 
National Forest (Appendix C). Prior to tlie 2006 
survey, tliis slug was known fi'om < 20 sites in 
Montana (Hendiicks 2003; Hendi'icks et al. 2006, 
2007) and remains documented at < 30 sites 
following tlie 2007 survey Canopy at tlie 2007 
sites included western redcedar, western hemlock, 
Engelmann spiiice, grand fir, Douglas-fir, alder, and 
black Cottonwood. Single animals were found at 
each site undei' downed rotten wood or rock. 

Smoky Taildiopp cr (Prophysaon humileY We 
found this species at six sites between 3140-4640 
ft elevation in Lincoln, Missoula, and Sanders 
counties, on the Libby and Thiee Rivers RDs of tlie 
Kootenai National Forest, and the Ninemile and 
Plains RDs of the Lolo National Forest (Appendix 
C). This species is known only from northern 
Idaho and adjacent nortliwestern Montana (Pilsbry 
1948: Frest and Johannes 1995, 2001; Hendricks 
2005; Hendiicks et al. 2006, 2007). Prior to 2004 
tliis slug was known in Montana from a single 
site. Following tlie 2007 surveys, it has now been 
documented at about 41 sites. Canopy at tlie 2007 
sites included western redcedar, giand fir, Douglas- 
fir, Engelmann spruce, western laixh, subalpine fir, 
lodgepole pine, alder, and black cottonwood. Up 
to thiee individuals (nine total) were found under 
downed rotten wood. 



Lyre Mantleslug (Udosarx lyraid)\ We found this 
species at one site at 4640 ft elevation in Missoula 
County, on the Ninemile RD of the Lolo National 
Forest {Appendix C). This species is restricted 
to nortliem Idaho and adjacent parts of western 
Montana (Webb 1959; Russell and Webb 1980; 
Frest and Johannes 1995, 2001; Hendricks 2003; 
Hendiicks et al. 2006, 2007). Two subspecies are 
described but we are unable to distinguish these, 
and assign our records only to tlie species level. 
Altliough known fioin Montana since 1965, theie 
remain only seven reported locations in tlie state, 
tliree of which were found in 2006. Globally, 
tliere are fewei' than 15 records (Hendiicks et al. 
2006, 2007). Canopy at tlie 2007 site included 
Engelmann spiiice, Douglas-fir, and alder. Only 
a single animal was found, under downed rotten 
wood. 

Sheathed Slug (Zacolens idot\oensis)\ We found 
tliis species at four sites between 2680-3630 ft 
elevation in Lincoln County, on tlie Libby, Rexford, 
and Tliree Rivers RDs of tlie Kootenai National 
Forest (Appendix C). This species is restricted to 
noitliern Idaho and adjacent nortliwestern Montana 
(Pilsbry 1948: Frest and Johannes 1995, 2001, 
Hendiicks 2003; Hendricks et al. 2006, 2007). The 
total number of documented Montana localities is 
33 following the 2007 suiTeys. Canopy at the 2007 
sites included western redcedar, Engelmann spruce, 
lodgepole pine, ponderosa pine, western larch, and 
aldei; One to two individuals were found at each 
site (six total), under wet downed rotten wood, rock 
fragments, or leaf litter. 



Conclusions and Recommendations 



Tlie number of new locations we discovered 
(1 16 total) for 1 6 of 1 9 species of land mollusks 
of conseiTation concern during our 2005-2007 
surveys in die Montana portion of die USPS 
Noitiiem Region aiea underscores our conclusion 
Uiat cunent knowledge of Uie distribution, ecology, 
and status of this suite of species remains far 
from thorough. Altliough we have accumulated 
many new records for Species of Concein on 
Forest Service lands during the suiTeys, many 
taxa remain documented from < 25 locations. Of 
course, some of these may be so raie diat they will 
never be found at more Ihan a handful of sites, 
but we predict new locations eventually will be 
documented for most. 

We diink, however, an additional yeai" of /io«- 
rmidom suiveys, similar to diose of 2005-2007, 
would help fill disti'ibution gaps and gather 
additional habitat infomiation. We also suggest 
a minimum of one additional yeai" of ra/nrhtn site 
surveys are needed to detemiine site occupancy 
rates as a measure of status in various habitats. It is 
important Uiese be conducted under environmental 
conditions when land mollusks are most likely to 
be active. 

Developing predictive habitat or ecological niche 
models may also prove usefiil for guiding surveys 
of some species groups, especially Uiose associated 
wiUi the moist forest types mentioned earlier. To 
increase (he utility of predictive habitat models, 
it is impoitant tiiat more-detailed habitat data are 
recorded when and where SOC and SOI species are 
found. Recent examples of the use of predictive 
models for conservation management of rare land 
mollusks in the Pacific Nortliwest and Black Hills 
are Dunk et al. (2004), Gaines et al. {2005), and 
Weaver etal. (2006). 

Appendix E provides an overview of inputs, an 
example output, and a guide to interpretation of 
outputs for the predicted disti'ibution models. 
Model output for individual species can be viewed 
under individual species accounts on die on-line 
Animal Field Guide website at http://fieldguide. 
iiit.gov . Overall, the predicted distiibution model 



output indicates that tlie known ranges of a number 
of species may be extended by focusing fiiture 
survey effort in the vicinity of Glacier National 
Pai'k, die Seeley/Swan aiea, die SouUi Fork of the 
Flatiiead River diainage, die Mission Mountains 
and die Bitterroot Mountains in Mineral County 
(Figme 2). Major gaps in survey effort on U.S. 
Forest Service lands still exist in a number 
of aieas, but model output indicates that die 
greatest number of gaps in known species ranges 
(especially Species of Concern) may be filled in 
witii additional suiveys focused on portions of 
noitliwestem Montana, the Bitterroot Mountains 
in Mineral and Ravalli Counties, Beaverhead, 
Madison, and Silver Bow Counties, and portions 
of the island mountain ranges in eastern Montana 
(Figure 2). In addition to diese gaps, a summaiy of 
elevations on U.S. Forest Service lands in Montana 
and elevations of mollusk sui^ey locations in 2005- 
2007 indicates diat high elevation sites have been 
under-sampled (Figure 3). We intend to target 
these areas during the 2008 field season. 

In addition to targeting diese ai'eas for survey, 
we continue to recommend die following actions 
as expressed in pievious reports (Hendiicks et 
al. 2006, 2007): (1) Survey and modeling efforts 
should continue to be coordinated with die Idaho 
CDC and MTNHP; this coordination is especially 
desired to detemiine more fidly die status of the 
many SOC and SOI species shared in the two 
states; (2) There remains a need for genetic studies 
to infonii cuiTent taxonomy of some species. We 
diink some taxa cuirently considered conspecific 
widi coastal populations (e.g.. Robust Lancetooth, 
Reticulate Taildropper, Blue-gray Taildiopper, and 
Papillose Taildiopper) may prove to be distinct 
sister species (see discussions in Leonard et al. 
2003, Ovaska et al. 2004), similai' to die results of 
recent genetic studies of some Pacific Noitliwest 
amphibian genera (G.g.,Asccq?hns\ Dicatuptodoii, 
Phthodoj]). Similai'ly, OreoheUx species across 
die region diat are similai^ in morphology may 
represent one species or numerous ciyptic species. 
A genetic study of die phylogeography of diese 
species is needed in ordei' to prevent loss of 
potentially unique species and/or populations. 



Finally, high elevation taxa like tlie Alpine 
Mountainsnail and Carinate Mountainsnail would 
benefit flom genetic analyses to detenniiie if tlieir 
isolated populations ai'e genetically identical, 
or if they ai'e distinct eiiougli to be classified as 
recently-isolated subspecies or sistei' species 
(see Appendix D): (3) Development of a field 
guide to tlie land mollusks of Montana would 
compile largely inaccessible pieces of information 



to heighten awareness of this overlooked and 
poorly understood group of animals, and provide 
biologists and managers some of tlie basic 
tools they need to make infomied management 
decisions. (4) Finally after development of a field 
guide, we tliink it would be usefiil to conduct 
workshops on land mollusk identification and 
management. 




Figure 2. PnonttzaUon (nmnbers m ovals) of areas lacking land mollusk surveys on U.S. Forest Service Northern Region lands 
as determined by predicted distribution models. Red ovals encompass areas most likely to residt in extensions of known species 
ranges. Green ovals encompass areas most likely to residt in the filling in of major gaps in known species ranges. 
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Figure 3. Distribution of elevations on U.S. Forest Service No/them Region lands in Montana across a 90-meter 
pixel digital elevation model (mean elevation indicated by dashed line). The distribution of land molhisk sui-vey 
locations is significantly lower than the distribution of available elevations as indicated by the siq?erimposed tall 
dark line (mean), short dark line (median), and 2 tall light lines (2 standard deviations above and below mean). 
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Appendix A. Global and State Rank Definitions 



Heritage Program Ravks 

Tlie international network of Natural Heritage Programs employs a standardized tanking system to denote 
global (range-wide) and state status. Species are assigned numeric ranks ranging from 1 to 5, reflecting 
the relative degiee to which they are ''at-risk''. Rank deflnitions aie given below. Anumber of factors are 
considered in assigning ranks — the number, size and disti'ibution of known ''occurrences" or popula- 
tions, population tiends (if known), habitat sensitivity, and threat. Factors in a species' life history Uiat 
make it especially vulnerable are also considered (e.g., dependence on a speciflc pollinator). 

Global Rank Definitions (NatureServe 2003) 
Gl Qitically impeiiled because of extreme rarity and/or other factors making it highly 

vulnerable to extinction 
G2 Imperiled because of rarity and/or othei' factors making it vulnerable to extinction 

G3 Vulnerable because of rai ity or restiicted range and/oi" otlier factors, even Uiough it may 

be abundant at some of its locations 
G4 Appaiently secure, though it may be quite rare in parts of its range, especially at the 

periphery 
G5 Deinonstiably secure, tliough it may be quite rare in parts of its range, especially at tlie 

periphery 
Tl-5 Infraspecific Taxon (tiinomial) — Tlie status of infraspecifictaxa (subspecies or 

varieties) are indicated by a "T-rank" following tlie species' global rank 

State Rank Definitions 

51 At high risk because of extremely limited and potentially declining numbers, 
extent and/or habitat, making it highly vulneiable to extirpation in tlie state 

52 At risk because of veiy limited and potentially declining numbers, extent and/or 
habitat, making it vulnerable to extirpation in the state 

53 Potentially at risk because of limited and potentially declining numbers, extent 
and/or habitat, even though it may be abundant in some areas 

54 Uncommon but not raie (although it may be raie in paits of its range), and usually 
widespread. Appaiently not vulnerable in most of its range, but possibly cause for 
long-tenn concern 

55 Common, widespread, and abundant (although it may be rare in parts of its 
range). Not vulneiable in most of its range 

Combination Ranks 

G#G# or S#S# Range Rank — A numeric range rank (e.g., G2G3) used to indicate uncertainty about 
Uie exact status of a taxon 

Qualifiers 

MR Not ranked 

Q Questionable taxonomy tliat may reduce conservation priority — Distinctiveness of 

tliis entity as a taxon at the current level is questionable; resolution of this uncertainty 
may result in change fiom a species to a subspecies or hybrid, or inclusion of diis taxon 
in anodier taxon, witli the resulting taxon having a lower-priority (numerically higher) 
conservation status rank 
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X presumed Extinct — Species believed to be extinct throughout its range. Not located 

despite intensive searches of historical sites and other appropriate habitat and virtually 
no likelihood tliat it will be rediscovered 

H Possibly Extinct — Species known from only historical occurrences, but may never-the- 

less still be extant: furtlier searching needed 

U Unrankable — Species cunently unrankable due to lack of infoimation or due to substan- 

tially conflicting infonnation about status or tiends 

HYB Hybrid — Entity not ranked because it represents an interspeciflc hybrid and not a species 

? Inexact Numeric Rank — Denotes inexact numeric rank 

C Captive or Cultivated Onfy — Species at present is extant only in captivity or cultiva 

tion, or as a reintioduced population not yet established 

A Accidental — Species is accidental or casual in Montana, in other words, infiequent and 

outside usual range, hicludes species (usually birds or butterflies) recorded once or only a 
few times at alocation. Afewof these species may have bred on the one or two occa- 
sions they were recorded 

Z Zero Occurrences — Species is present but lacking practical conservation concern in 

Montana because there are no definable occurrences, altiiough the taxon is native and 
appears regulaily in Montana 

P Potential — Potential that species occurs in Montana but no extant or historic occurrences 

are accepted 

R Reported — Species reported in Montana but without a basis for eiUier accepting or 

rejecting the report, or Uie report not yet reviewed locally. Some of Uiese aie very recent 
discoveries for which Uie progiam has not yet received first-hand information; others aie 
old, obscure reports 

SYN Synonym — Species reported as occurring in Montana, but the Montana Natural Heritage 

Progiam does not recognize the taxon; therefore the species is not assigned a rank 

* A rank has been assigned and is under review. Contact tiie Montana Natural Heritage 

Progiam for assigned rank 

B Breeding — Rank refers to tiie breeding population of the species in Montana 

N Nonbreeding — Rank refers to Uie non-breeding population of the species in Montana 
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Appendix B. Data Forms 



Data Form for Terrestrial Mollusk Site Surveys 

Locality Inform a lion 



Ecoregion 


Sample Black 


Sits Ho 


Locality 


Slate 


County 


Map 
Name 


T 


R 


S 


Sectian 
Description: 


Owner 


Map 
Elevation 


FT 


Datum. 


UTM 
Zone 


UTM 
Ea^t 


UTM 
Worth 



Habitat Information 



Dale: 



Ob*;erver(:o) 



Begm 
Time 



End 
Time 



Total Person 
Minules of Search 



Area (M') 
Searched 



Percentage of Site Searched 
1-25 26-^U 51-75 7fi-IOO 



Percenl 
Slope: 



Asp eel N NE NW 



SE SW 



W 



Habilat Type 

Spring/Seep Slreamade Talus Deciduous Forest Conifer Forest Mixed Forest Shiub/Sleppe Grassland 
Other 



Primary Canopy Species 



Overall Percent Canopy Cover 

1-25 26-SU 51-75 76-100 



Canopy Species Average DBH (cm) 
0-5 5-15 15-30 3U-6U >60 



PholD Frame Number (s) 
/ Desi;nption(:o) 



Weather Clear Partly Cloudy Overcast Rain 



Snow 



Air Temp 



Soil Temp: 



Soil Moisture 

Dry Damp Wei Standing Water Snow 



Rock Type" Igneous Metamorphic 
Note Specific Type (e g. limestone, granite) 



Sedimentary 



Habilat 
Threats 



MoUusk Species Information 



Species 



Number Alive and/or Dead, Size, and Time at First Detection (e g, 2 alive &4 deads 1 5mm Diameler or TL@ 10 minutes) 



Tissue Number (eg, HOOIA) 



Voucher Number 
Sl Description 



Substrate Assoaation (Circle) 
underwood under 4 -2D cm rock fragments under >2Dcm rock fragmenls 
under biy op byte mal on bryophyte mat in rock fracture Other 



Species 



Number Alive and/or Dead^ Size, and Time at First Detection (e g, 2 alive & 4 deadx 1 5mm Diameter or TL@ 10 minutes) 



Tissue Number (eg, HOOIA) 



Voucher Number 
Sl Description 



Substrale Association (Circle): 
underwood under 4 -2 0cm rock fragments under >20cm rock fragments 
under bryophyte mat on biynphyte mat m rock fracture Other 



Species 



NimiberAhve and/or Dead, Size, and Time at First Detection (e g, 2ahve&4 deadx 1 5mm Diameler or TL@ 10 minutes) 



Tissue Number (e g, LCQOIA) 



Voucher Number 
& Description 



Substrate Association (Circle): 
underwood under 4 -2 0cm rock fragments under >20cm rock fragments 
under bryophyle mat on bryophyte mat in rock fracture Other 



Species 



Number Ahve and/or Dead, Size, and Time at First Detection (e g, 2ahve&4 deadx 1 5mm Diameler or TL@ 10 minutes) 



Tissue Number (eg, GOOIA) 



Voucher Number 
life Description 



Substrale Association (Circle)" 
underwood under 4 -2 0cm rock fragments under >20cm rock fragments 
under bryophyle mat on bryophyte mat in rock fracture Other 
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Site Map for Terrestrial Mollusk Site Surveys 



Grid 


Scalr: 
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* Draw a rough sketch of the site labehng major features such as sticams, talus slopes, habilat cover types, etc. 
Be suie to mdicate where amuials were detected and label the foUowmg locations ou the map: G = GPS reading, 
and P-^ = photo locations and directions of photos. 
Other Notes: 
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Site ID (ecoregioii, sainplr block, site number) 



Date: 



Mollnsk Species Information Continued 



Species 



Niimber Alive and/oi" Dead, Size, and Time alFirslDelectian (eg , 2 alive &4 deads 15mm Diameter ar TL@ 10 minutes) 



Tissue Number (e g, HOOIA) 



Voucher Numb er 
&r Descnption 



Substrate Associalian (Circle) 
underwood under 4-20cm rock fragments under >20cm rack fragments 
under biyaphyte mat an bryophyte mat in rock fracture Otlier 



Species. 



Number Alive and/or Dead, Size, andTimeatFiislDelectian(eg , 2 alive &4 daids 15mmDiamelei"ar TL@ 10 minutes) 



Tissue Number (e g, HOOIA) 



Voucher Number 
&r Description 



Subsitrate Association (Circle) 
underwood under 4 -2 0cm rack fragments; under >I10cmrock fragments 
under bryophyte mat on bryophyte mat in rack fracture Other 



Species: 



Number Ahveand/or Dead, Size, and Time at First Delectian (eg , 2 alive &4 daids 15mm Diameter ar TL@ 10 minutes) 



Tissue Number (eg, LCOOIA) 



Voucher Number 
& Description 



Substrate Association (Circle) 
underwood under 4 -2 0cm rack fragments under >?Ocmrock fragments 
under bryophyte mat on bryophyte mat m rack fracture Other 



Species: 



Number Ahveand/or Dead, Size, and Time at Fust Detectian (eg , 2 alive &4 daids 15mm Diameter ar TL@ 10 minutes) 



Tissue Number (e g, GOOIA) 



Voucher Number 
Sl Deiocnption 



Substrate Association (Circle) 
underwood under 4 -2 0cm rack fragments under >30cmrock fragments 
under bryophyte mat on bryophyte mat in rack fracture Other 



Other Notes 



Otlier Sp eeies Information 



Other Species: 

(millipedes etc ) 


Time at First 
Detection 


Voucher 
Number 


Voucher Descnptian / Comments 


Other Species: 

(millipedes etc ) 


Time at First 
Detection 


Voucher 
Number 


Voucher Descnptian / Comments 


Other Species: 

(millipedes etc ) 


Time at First 
Detection 


Voucher 
Number 


Voucher Descnptian / Comments 


Other Species: 

(millipedes etc ) 


Time at First 
Detection 


Voucher 
Number 


Voucher Descnptian / Comments 


Other Species: 

(millipedes etc ) 


Time at First 
Detection 


Voucher 
Number 


Voucher Descnplian / Comments 


Other Species: 

(millipedes etc ) 


Time at First 
Detection 


Voucher 
Number 


Voucher Descnptian / Comments 
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Site Inform arion 

Ecoregion: One of the 14 ecoregion liections in Monlana or 6 in the Idaho Panhandle 

Sample Block: Identify three digit number of the sampling block, {range 001-999). 

Site No: Identify three digit number of the liite being surveyed within each sampling block (range 001-999) 

LocaUty: Describe the specific geographic location of the site so that the type of sile is described and the straight-hne air distance 

from one or more permanent features on a 7 5-minute (1 24,000 scale) topographic map records the po^ilion of the sile (e g , Large 

talus slope 1 5 miles north of Engl e Peak.N side of FS Road 225) 

State: Use the two-lelter abbreviation 

County: Use the full county name. 

Map Name: List the name of the USG5 1 5-minule (1 24^000 scale) topographic quadrangle map 

T: E.ecord the Township number and whether il is north or south 

R: E.ecord the E.ange number and whether it is east or west 

S: Record the Section number 

Section Description: Descnbe location of the site at the Vd of ^^ section level (e.g , SENE indicates SE comer of NE corner). 

Ownei-; Use abbreviation of the government agency responsible for managing the land you surveyed, (e.g USFi^, BLM) If private 

land was surveyed list the owner's full name to indicate that you did not trespass 

Map Elevation: The elevation of the site as indicated by the topographic map in feet (avoid using elevations from a GPS) 

Datnm: The map datum used (typically NAD 27 if off topographic map or WGSS4 if off GP;5 unit on standard setting) 

UTMZone: Universal Transverse Merc ator zone recorded on the topographic map 

UTMEa^: Universal Transverse Mercator easting coordinate in meters as recorded on the topographic map or GPS receiver Be 

sure to note any major differences between UTM coordinates on the map and those on the GPS receiver 

UTMNorth: Universal Transverse Mercator northing coordinate in meters as recorded on the topographic map or GPS receiver Be 

sure to note any major differences between UTM coordinates on the map and those on the GPS receiver 

Survey Inform ation 

Date: Use mi-DD-YY format (e g 05/12/00 for May 12 of 2000). 

Observers: List names or initials of individuals involved with survey of this site and circle the name of the recorder 

Begin Time: List the time the survey began in 24-hour format 

End Time: List the time the survey ended in 24-hour format 

Total Pel' son Minutes of Search: B.ecord the total person minutes the site was searched (e g if one person surveys for 15 minutes 

and another suiveys for 30 minutes, but takes 5 minutes to measure a specimen the total person minutes is 40 minutes) 

Area (M ) Searched: Aiea in square meters that was surveyed 

Percent of Site Searched: Circle the appropriate category 

Percent Slope: Percent slope of site Enter range if variable 

Aspect; Circle piimaiy aspect of the site. 

Habitat Type: Circle the appropriate habitat type. 

Primary Canopy Species: List the major plant species in the canopy (e.g , red cedar, western hemlock, grand fir, ninebark) 

Overall Percent Canopy Cover; Circle the appropriate category for total canopy cover 

Canopy Species Average DBH: Circle the appropriate category 

Photo Frame Numb er(s) / Descriptions: The number of the photo as viewed on the camera^ s view screen and a description of the 

contents of the photograph (e.g , #13 = 1 x Oreohehy stiigosa and #14 -IS = 5 s habitat) Take photos of all portions of the site and 

anything else that may be of interest (e g , millipedes, potential site threats) 

Weathei" Circle weather condition during survey 

Au- Temp: Record air temperature in °C at chest height in the shade °C = (°F-32)/l 8 

Soil Temp: Record soil temperaUire in °C at 10 cm depth °C = (°F - 32)/l 8 

Soil Moisture: Circle the appropriate category 

Rock Type; Circle the appropriate category, note specific type if known. 

Habitat Threats: Note impacts from grazing, logging, mining, flooding, road building, weeds, tire, etc. 

Sp ecies Inform ation 

For each species, record the genus name and species, if known If species cannot be identified in the field, place abnef description of 

their morphology here Record the number alive and dead, and size range for individuals encountered, and time at first detection for 

the first individual encountered (e g , 2 x 15 mm diameter (shells) or TL = S0-90mm (slugs) @ 10 minutes) Record the tissue 

number or range of tissue numbers for tissue samples collected (see tissue collection protocols) Record the preliminary museum 

voucher specimen number and description for voucher specimens collected (see voucher specimen collection protocols) Circle the 

substrate the animal was associated with at time of detection Record the presence of other species detected at the site (e g , 

milhpedes), the time at first detection, and the voucher number and description of animals collected (see voucher and tissue collection 

protocols) 



Appendix B - 4 



Appendix C. USFS Northern Region Survey Sites in 2007 

FOR Land Mollusks 



O 



Forest^ 


County 


NAD 27 


Site Name 


Elev 
(ft) 


Date 


SOC/SOITaxa"* 


KO 


Lincoln 


N48.92831 W114. 74961 


Lewis Creek 


4500 


6 Sep 


Ma my 


KO 


Lincoln 


N48.8911S W114.79566 


Clarence Creek 


3940 


6 Sep 


Mamy 


KO 


Lincoln 


N48. 86820 W114. 31010 


Cat Creek 


3600 


6 Sep 


^ 


KO 


Lincoln 


N48. 81930 W114. 88702 


Grave Creek CG 


3040 


7 Sep 


:+ 


KO 


Flathead 


N48.63106W114. 74650 


Lower Sunday Creek 


3400 


7 Sep 


* 


KO 


Flathead 


N48.54058W114. 86129 


Upper Sunday Creek 


4810 


7 Sep 


:+ 


KO 


Flathead 


N48.54771W114. 89754 


Upper Basin Creek 


4640 


7 Sep 


^ 


KO 


Lincoln 


N48.59669 W115.67039 


Loon Lake CG 


3630 


1 Oct 


Kobii,Prhu, Zaid 


KO 


Lincoln 


N48.95140W115. 88522 


West Fork Yaak River 


4605 


2 Oct 


Heca, Kobii 


KO 


Lincoln 


N48.91225W115. 86121 


Beetle Creek 


4800 


2 Oct 


Heca, Kobu 


KO 


Lincoln 


N48.93122W115. 72604 


West Fork Yaak River Falls 


3360 


3 Oct 


Kobii, Zaid 


KO 


Lincoln 


N48.92445W115. 53957 


Solo Joe Creek 


4020 


3 Oct 


Heca, Kobii.Prliu 


KO 


Lincoln 


N48.75180W115. 37236 


Little North Fork Big Creek 


2680 


3 Oct 


Kobii, Zaid 


KO 


Lincoln 


N48. 87092 W115. 25772 


Lowei Sullivan Creek 


2980 


4 Oct 


Heca, Kobii, Zaid 


KO 


Lincoln 


N48.92896W115. 23691 


Poverty Creek ponds 


3010 


4 Oct 


:+ 


KO 


Lincoln 


N48.29531 W115.6313S 


Granite Creek trailhead 


3140 


5 Oct 


Kobu, Many, Prhii, Raab 


LC 


Lewis & Clark 


N47 .33039 W112. 73004 


Welcome Guard Station 


5920 


22 Jul 


* 


LC 


Lewis & Clark 


N47 .32235 W112. 75405 


Upper Dearborn River 


6110 


22 Jul 


:+ 


LC 


Lewis & Clark 


N47 .33156 W112. 82243 


Upper Halfinoon Creek 


6950 


22 Jul 


:+ 


LC 


Lewis & Clark 


N47 .33390 W112. 82540 


^^Halfinoon Ridge" 


7150 


23 Jul 


^ 


LC 


Lewis & Clark 


N47 .33552 W112. 84132 


Green Fork 


6570 


23 Jul 


^ 


LC 


Lewis & Clark 


N47 .33130 W112. 84 368 


Green Fork wetland 


<j670 


23 Jul 


^ 


LC 


Lewis & Clark 


N47 .31065 W112. 85110 


Scapegoat Mountain west 


8075 


24 Jul 


Oral 


LC 


Lewis & Clark 


N47 .30867 W112. 85936 


Scapegoat Plateau 


7960 


24 Jul 


Orel 


LO 


Missoula 


N47. 18370 W113. 51700 


Camp Paxson 


4000 


20 Sep 


* 


LO 


Missoula 


N46.74104W114. 54756 


^^NorthFork Granite Creek'' 


4400 


24 Sep 


Raab 


LO 


Miiieial 


N46.78493 W114.62215 


"North Fork Oriole Creek" 


4960 


24 Sep 


Raab 


LO 


Mineral 


N46. 83879 W114. 71236 


Suiveyors Creek 


3900 


24 Sep 


Raab 


LO 


Mineral 


N46.90919W114. 80193 


Clearwater Crossing CG 


3470 


24 Sep 


^ 


LO 


Missoula 


N47. 13948 W114. 55600 


Monciire Creek 


3960 


11 Oct 


* 


LO 


Missoula 


N47 .21630 W114. 58950 


Camp Creek 


4280 


11 Oct 


Prliu, Raab 


LO 


Sanders 


N47 .27145 W114. 72590 


Upper Siegel Creek 


3920 


11 Oct 


Prliu, Raab 



Forest^ 


County 


NAD 27 


Site Name 


Elev 


Date 


SOC/SOITaxa"* 


LO 


Misfe'oula 


N47 .25455 W114. 73360 


Elista die Creek 


4640 


11 Oct 


Prim, Utlly 


LO 


Miueial 


N47 .28139 W115. 34135 


Upper Ward Creek 


4100 


13 Oct 


Heda, Kobu, Maniy, Popo 



" National Forests KO (Kootenai), LC (Lewis and Clark), LO (Lolo) 

* = only non 30C/30I moUiisks detected; 30C/S0I taza codes Heca (Pale Jumping-slug, Hemphdha camslns). Heda (Marbled Jump ing-s lug, Hemphdha danislsi), Kobu 
(Pygmy Slug, Kootenai bnrkei), Mamy (Magnum Mantle-slug, Magmpelta mycophaga). Oral (Alpine Mountainsnail, Oreohehx aipind), Orel (Carinate Mountainsnail, Oreoheltx 
elrodi), Popo (Humped Coin, Polygyrella polyg}'!-eUa)^ Prhu (Smoky Taildropper, Prophysaon humile)^ Raab (Fir Pinwheel, Radiodisciis abietum), Udly (Lyre Mantleslug, Udo- 
sarx lyratd), Zaid (Sheathed Slug, Zacoleus idahoensis 



O 



K> 



Appendix D. Observations of the Scapegoat Mountaevsnails 

{Oreohelix) 



OBSERVATIONS OF THE SCAPEGOAT MOUNTAINS NAILS (OREOHELIX) 

During an ascent of Scapegoat Mountain (9202 ft elevation) on 24 July 2007 (Rocky Mountain RD, 
Lewis and Clark National Forest) byHendi'icks and companions, and explorations of the Scapegoat 
Plateau about I.0-L5 miles west of tlie summit (7800-8100 ft elevation), two species ofmountainsnail 
(OreohelLx) were found in bairen limestone alpine teirain in tlie basin (cirque) di'aining noitli into Green 
Fork creek. Collections of shells and live animals were made in tlie field for later determination, but the 
two species appeared to represent tlie Alpine Mountainsnail {OreoheUx alpifta) and tlie Caiinate Moiin- 
tainsnail (Oreohelix elrodi), or an undescribed fonn closely related to OreoheUx elrodi. 



Suiprisingly of the 12 sites occupied by snails and searched for at least 20 minutes, O. alpi}ia?a\A O. 
elrodi were present at half each (shells and live individuals at all sites) but both species were never found 
togetlier, nor did the ai'easofoccunence for both species overlap (i.e., appaiently no aiea of sympatiy). 
Tliis result is statistically significant (Fisher Exact Test, P = 0.0137) assuming equal probability of finding 
one and two species at an occupied site. O. cdpitm was found in talus sites under loose rocks about 1 .0- 
1.25 miles west of the summit of Scapegoat Mountain (Figures D1-D4), whereas O. ^//W/ was present 
in talus and more stablerocky sites under rocks 1.25-1.5 miles west of tlie summit (Figures D5-D8). The 
full extent of occurrence for botii species on the Scapegoat Plateau remains unknown, because only the 
central region of tlie plateau was searched and investigation of this area was very incomplete. 




Figure DL Alpine Mountainsnail (Oreohe'ix aipwa) on 24 
July 2007 in situ under limestone block at about 8000 fi on 
the Scapegoat Plateau, 




Figure D3. Local habitat typical of sites occupied by the 
Alpine Moimtainsnail on the Scapegoat Plateau. 




Figure D2. Alpine Moimtainsnad on the underside of a 
limestone block. 
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Figure D4. Landscape view from 8000 ft of the Scapegoat 
Plateau where Alpine Moimtainsnails were found. Summit 
of Scapegoat Mountain (9202 ft} in the background. 
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Figure D5. Carmate Mountainsnail (Oreohelix elrodi) on 
24 Jidy 2007 in situ under limestone block at about 7900ft 
on the Scapegoat Plateau. 




Figure D6. Carmate MowT.tainsnails. 
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Figure D7. Local habitat typical of sites occiq>ied by the 
Carinate Moiintainsnad on the Scapegoat Plateau. 



■- '^±'. 



Figure D3. Landjc:pe vv,.-.. ;^j^^/ SOOO ft cftih^ Sccpegcct 
Plateau where Carmate mountainsnails were found. Occu- 
pied sites were along the margins of the relatively flat basin 
in the center-right of the photograph. 



Alpine Mountainsnail (Oreohelix alpind) 

Preliminaiy analysis ofmoiphological data show shell diameter and height of the Scapegoat Oreohelix 
alpina are consistent with shell morphology data from the otlier known locations for tlie Alpine Moun- 
tainsnail in the Mission and Swan ranges to the west (Table Dl: see also EIrod 1903b; Hendricks 1998; 
Fairbanks 2002). Mean values for shell tiaits of the Swan Range individuals are laiger dian for the Scape- 
goat and Mission specimens because of a lack of immature animals in the Swan Range sample. 

Table Dl. Shell morphology for Oreohelix cdpinafrom three Montana moimtam ranges. HD is the ratio of height divided by 
diameter. Values are means (SD); range. 



Variable 


Scapegoat Plateau (62) 


MissioD Mountains (17) 


Swan Range (10) 


Diiun (mm) 


7.7(1.4); 2.4-11.7 


7.4 (2.1); 2.6-9.6 


9.7(1.0); S. 5-11.3 


Heiglit {mm) 


4.6(0.8); 14-7.0 


4.1 (1.3); 1.3-5.2 


5.4 (0.5); 4.9-6.4 


HD 


0.6 (0.04); 0.5-0.7 


0.5 (0.03); 0.5-0.6 


0.6 (0.3); 0.5-0.6 


No. whorls 


3.9(0.5); 1.6-4.6 


3.6 (0.6); 2.0-4.1 


4.1 (0.2); 3.8-4.3 



Plotting shell height against shell diametei" sliows a relatively linear relationship for the combined sample 
from all three mountain ranges (Figure D9), witli no obvious areas of segregation in the samples. 
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Oreohelix alpina 
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Figure D9. Relationship of shell height to shell diameter for the Alp me 
Monntamsnail (Oreohelix alpinaj. Sample includes specimens from the 
Scapegoat Plateau (62), Mission Mountains (17), Swan Range (10). 

Carioate MouDtaiosnail (OreoheOx elrodt) 

Preliminaiy analysis of morphological data show shell diameter and height of the Scapegoat Oreohelix 
elrodt are much sinallei" than shells fiom the otlier known locations foi' tlie Carinate Mountainsnail in the 
Mission and Swan ranges to the west (Table D2; see also Ehod 1903a; Fairbanks 1984; Hendricks 1998). 
Mean values for shell traits of individuals from die Mission Mountains and Swan Range ai'e larger than 
for tlie Scapegoat, in part because of a lack of inunature animals in the Missions saiuple, but also because 
no animals found on the Scapegoat Plateau exceeded 15 mm diameter' whereas 27 of 36 from tlie com- 
bined Mission/Swan sample exceeded tliis value. Lack of laiger individuals in tlie Scapegoat sample does 
not appeal' to be an artifact of sampling. Over 500 shells on the Scapegoat Plateau were encountered, 
and die laigest individuals were collected fr'oni each site where tliey were present. Presence of few small 
animals in the combined Mission/Swan sample, however, indicates under-sampling of diat size class. 

Table D2. Shell morphology for Oreohelix elrodifi-om three Montana mountain ranges. HD is the ratio of height divided by 
diameter. Values are means (SDj; range. 



Variable 


Scapegoat Plateau (51) 


Mission Mountains (20) 


Swan Range (16) 


Di^n (mm) 


10.5 (2.5); 3.2-14.7 


18.0 (9.7); 9.7-22.7 


14.4(7.8); 3.1-23.7 


Heiglkt (mm) 


4.7(1.1); 1.5-6.4 


8.7 (1.7); 4.5-11.4 


7.1 (4.1); 1.4-12.7 


HD 


0.5 (0.3); 0.4-0.5 


0.5 (0.03): 0.4-0.6 


0.5 (0.04); 0.4-0.6 


No. whorls 


3.7(0.5); 1.8-4.5 


4.7(0.4); 3.7-5.3 


4.1 (1.3); 2.1-5.5 



There are at least two possibilities explaining why die Scapegoat OreoheUx elrodt are physically smaller 
tlian individuals from other populations. First, the alpine habitat occupied by this species on the Scape- 
goat Plateau is unlike diat of the otlier populations in tlie Mission and Swan ranges, where occupied talus 
slopes are well below treeline and imbedded in amatiix of tiees andshmbs (Hendiicks 1998; Hendi'icks 
2003). Fuithennore, the rocky substrate is different; the Scapegoat population was found amongst linie- 
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stone talus, while the species is found in diorite or argillite talus in the Missions and Swans. TIius, a more 
favorable climate and habitat in tlie Missions and Swans could favor a largei' average adult body size. 

Second, the Scapegoat Oreobelix is not O. elmdi but a closely related undescribed tax on. As mentioned 
above, it's presence above tieeline in limestone terrain is unlike tliat for other known populations of O. 
elrorlt. At the very least, tlie Scapegoat population could represent a distinct subspecies ofO. elrodi, and 
tlie description of such could be justified on the basis of shell morphology alone. 

We prefer a conservative approach until genetic and anatomical studies can be completed on tlie Scape- 
goat population and contrasted witli the Mission Mountains and Swan Range populations. Tiie relation- 
ship of shell height to shell diameter of a sample using individuals from all three mountain ranges is 
lineai\ suggesting the Scapegoat population is on the same growth trajectory as tlie otlier populations 
(Figure DIO). Tlius, aldiougli gene flow is no longer occuiring between these populations, Uiey could still 
be moiphologically and genetically indistinguishable, as is the case with widespread populations of Oreo- 
lielix cooperi in eastern Montana and tlie Black Hills of Wyoming and Soutii Dakota {Weaver et al. 2006; 
Anderson and Schmidt 2007; Anderson et al. 2007). 





Oreohelix eirodi 
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Figure DIO. Relationship of shell height to shell diameter for the Carinate 
Monntamsnad (Oreohelix elrodi}. Sample includes specimens from the Scapegoat 
Plateau (51), Swan Range (16), and Mission Mountains (20). 
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Appendix E. Predicted Distribution Models for Northern 

Region Land Mollusks 



Introduction 

Tliis appendix gives background infomiation on the modeling effort. Model output for all terresti'ial mol- 
lusks is posted under tlie individual species accounts on the online Montana Field Guide at: http://fleld- 
guide.mt.gov/ 

Despite recent survey efforts during the 2005 thiough 2007 field seasons, observation records for tenes- 
trial mollusks in Montana are very limited for most species {Table Al). Duetotlie lackof obsei'vation 
records, predicted disti'ibution models were created in order to: (1) better understand the potential ranges 
of species in Montana; (2) identify areas of predicted distribution tliat should be targeted for ftiture sur- 
veys; (3) detennine which enviionmental variables have the highest predictive association with presence 
of individual species; and (4) identify spatial patterns in the suitability of habitat for tliese species. 

Modeling Approach 

Program Maxent Version 2.3 (Jaynes 1957; Miroslav et al, 2004, 2005, 2007; Phillips et al. 2004, 2006; 
www.c5.princeton.edu/-schapire/maxent/) was used to create statewide predicted disti'ibution models 
based on locations of observations for individual species in tlie statewide Point Observation Database and 
6 continuous {elevation, slope, landsurface curvature, annual precipitation, maximum July temperature, 
minimum January temperature) and 5 categorical {aspect, geology, 1 992 National Land Cover Data, soil 
temperature class, STATSGO soils) enviionmental layers. Observation locations for each species were in- 
put in comma delimited format and environmental layers were fitted to 90 meter grid cell pixels in ASCII 
format. Botli species obseivation locations and enviionmental layers were projected to Nortli American 
Datum 1983 in Montana State Plane. 

Program Maxent Version 2.3 (Jaynes 1957; Miroslav et al. 2004, 2005, 2007; Phillips et al. 2004, 2006; 
www.cs.princeton.edu/-schapire/maxent/) generates probability distributions for all enviionmental layers 
over all pixels that are associated witli observations of tlie species. Tliese environmental ''features" are 
used to constiain tlie probability distribution computed to predict the statewide distribution of species so 
tliat estimated distributions match chaiacteristics of the empiiical distributions resulting flom the positive 
data. Depending on the amount of positive data available the estimated distiibution is automatically con- 
strained to match tlie empirical average (linear feature), average and variance (quadjatic feature), covai'i- 
ance (product feature), proportional occunence (tlireshold feature), or average below and constant above 
a certain point (liinge feature). Categoiical environmental variables have a discrete feature foi' eveiy 
possible value of the variable so that tlie estimated disti'ibutions have the same proportional representa- 
tion of each categorical value. In practice Maxent 2.3 avoids oveifitting of models to the training data by 
^'regulaiizing" or relaxing the feature constraints so that feature expectations only have to be close to tlie 
empirical distributions ratlier tlian exactly equal to tliem. Aregularization parameter can be enteied in 
Maxent 2.3 in order to delineate just how close of a fit is needed between empirical features and estimated 
distributions. 

Maxent makes use of tlie Gibbs distribution which takes the fomi: 

P(x) = exp(cl *fl(x)-^c2*f2(x) + c3*f3(x)...)/Z 

Where cl, c2, c3, ... aie weighted constants, fl,f2, f3, ... aie tlie consti'ained enviromnental features, and 
Z is a scaling constant tliat ensures tliat P sums to 1 over all grid cells. Tlie Maxent 2.3 algoridim is guar- 
anteed to converge to values of cl, c2, c3, ... that give tlie unique optimum distiibution P, and, tlierefore, 
tlie outputs are deterministic. Tlie program begins with a unifomi distiibution (tlie distiibution witli maxi- 
mum entropy) and successively modifies each weighted constant on each iteration until either tlie change 
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in the ''gain" (Uie log of tlie number of grid cells minus tlie average of the negative log probabihtiesof die 
sample locations) falls below a settlireslioldor a set maximum number of iterations ai'e peif omied. Tlius, 
tlie Maxent 2.3 algorithm finds the most unifomi distribution (the one witli tlie maximum enti'opy) subject 
to the constraints of the environmental features. Maxent 2.3 can use eveiy giid cell tliat has values for 
all the environmental variables to calculate the disti'ibution. However, because there are a lai'ge number 
of 90-ineter giid cells for Montana (10,204 columns and 5,892 rows), and tlie modeling peifomiance of 
Maxent 2.3 does not improve significantly with very laige numbers of pixels, 60,000 backgiound pixels 
(added to tliose witli positive species data) were used to represent tlie variety of enviionmental conditions 
present in the data. Tliis limitation also serves to generalize or ^'regulai'ize" tlie model to make sure tliat it 
is not overfit to tlie limited number of locations where species were obsei^ed. 

Tlie main output is in the form of a cumulative fimction with grid cells having values of 1 00 representing 
tlie most suitable habitat (best fit to the positive data) while tliose witli values closer to are less suitable 
habitat (worst fit to the positive data). Note that Maxent Version 2.3 output is not inteipreted as prob- 
ability of occupancy. Other standard model output includes: (1) an evaluation of omission enor rates 
for tiaining and test data as a function of the cumulative tlireshold and overall predicted aiea. Note that 
all data available for teirestiial inollusks was used to train the models so tliere was insufficient data to 
independently test tlie models; (2) a Receiver Operating Chai'acteiistic (ROC) curve which evaluates the 
overall predictive power of tlie model; (3) a table of potential thiesholds which could be used to catego- 
rize the continuous output into a binomial predicted or not predicted output; (4) a statewide map showing 
the continuous probability function as well as the training and test data points; (5) response curves for 
individual enviionmental variables showing how each affects tlie output while all other environmental 
vaiiables aie held constant at their average sample value; and (6) jackknife charts showing the relative 
importance of all enviionmental variables as a function of tlie fit of tlie model to the data witli their exclu- 
sion or sole inclusion in the model. 

Below we provide an overview of terrestrial mollusk observations used in tlie modeling effort as well 
as comments on tlie model output. We then describe tlie enviromnental input layers and go tluough an 
example inteipretation of model outputs. Complete model output for individual species has been posted 
to individual species accounts on tlie online Montana Field Guide at: http://fieldgiiide.mt.gov/ 
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Table L Ovei-view ofTeirsstrial Molhisk Observations Used in Modeling Effort and Model Comment. Of the 58 terrestrial molliisks documented in Montana, 39 species had suf- 
ficient numbers ofobsei-vation records for modeling their predicted distribution and of these, II had poor fitting models because they were based on very few records. The table 
below siiinmarizes observation records for these species and provides a brief commentaiy on model limitations and model output. In general models performed well in terms of 
omission and commission error rates where positive data is present in a region. However, there is reason to believe that the models mayperfonn less well in regions that lack sur- 
veys (i.e., they are probably over fit to regions that have positive data from suiveys). All model output should be regarded as a first iteration with additional modeling needed (eg., 
the exact characteristics of the soils and geology classes need to be identified when these variables are found to be driving the models) 



t 



S cientif ic Name 

(Bold = ^Fcit:; of Concern) 


GRank/ 
SRank 


Total 
No. 

Records 


No. Records 

< 400 Meters 

Locational 

Uncertainty 


No. Spatially 

Unique Records 

Used to Train 

Model 


Model Comments 


SNAILS 1 


Allo^ona ptychophora 


G5/35 


42 


41 


35 


Model output indicateii surveys are most needed in and ai'ound Glacier National Park 
and Bitteiroot di^ainai^es of Ravalli County 


Angiiispira kochi 


G5/35 


87 


87 


74 


Model output indicateii surveys are most needed in and ai'ound Glacier National Park, 
the 3outh Foi'k of the Flathead River drainage and the Seeley/3wan ai'ea. 


Catmella rehdsri 


G3Q/SNR 


"\ 


T 


T 


Model based on exti'emely small number of observation records Model output 
indicates surveys aie most needed in the Big and Little Snowy, Judith, and Piyor 
Mountains. 


Cochlicopa lubrica 


G5/SNR 


1 


1 





Not modeled due to limited numbei' of observation records 


Columella columella 


G5/3NR 


■h 


1 


-1 


Model based on exti'emely small numbei' of observation records 3pecies likely went 
undetected a number of suivey sites due to its small size and low detection probability 


Columella edentida 


G5/3NR 


4 


4 


4 


Model based on exti'emely small numbei' of observation records 3pecies likely went 
undetected a number of suivey sites due to its small size and low detection probability 


Ciyptomasttc midlani 


G4/34 


65 


64 


58 


Model output indicates surveys are most needed in and ai'ound Glacier National Park, 
the 3outh Foi^ of tlie Flathead River di'ainage, the Seeley/Swan area, and the 
Bitteiroot di'ainages of Ravalli County. 


D'iscus bmnsoni 


Gl/Sl 


9 


4 


1 


Not modeled due to limited numbei' of observation records from same location 


Discus shimekii 


G5/S1 


10 


6 


6 


Model based on exti'emely small number of observation records Model output 
indicates surveys are most needed in southwest Montana 


Disciis whitneyi 


G5/35 


113 


113 


110 


Species IS clearly widespread 


Eucomdi^ fidvus 


G5/35 


103 


108 


107 


Species IS clearly widespread 


Hapiotrema vancouverense 


G5/S1S2 


]] 


10 


10 


Model output indicates surveys are most needed in Pi'ospector Creek west of 
Thompson Falls and Bitterroot drainages in Mineral County 


Havvaua minuscide 


G5/SNR 


1 


1 





Not modeled due to limited numbei' of observation records 


Microp}^sida ingersolli 


G5/35 


88 


88 


81 


Species IS clearly widespread Additional surveys in southwest Montana would seem 
likely to detect their presence 


Nesovitrea binneyana 


G5/35 


41 


41 


37 


Species is clearly widespread 


Nesovitrea electrina 


G5/SNR 


1 


1 





Not modeled due to limited numbei' of observation records 


Oreohelix atpina 


Gl/Sl 


9 


7 


9 


tlodel output indicates surveys are most needed in and ai'ound Glacier National Park, 
Mountains in the Bob Marshall Wilderness, and the Swan and Mission Mountain 

Ranges 



t 



Oreohelix amariradix 


G1G2/S1S2 


•<< 


5 


3 


Model based on extremely small numbei' of obseivati on records 3ur/eys conducted 
in tlie Lolo Creek di'ainage are most likely to detect the species in new localitieii. 


Oreohelix carinifera 


Gl/Sl 


12 


"\ 


3 


Model based on extremely liinall number of obijervation recordii Model output 
indicates sur/eys are needed on a vanety of BLM andU3F3 lands in drier habitats 
that have not yet been targeted for surveys 


Oreohelix elrodi 


Gl/Sl 


28 


6 


6 


Model based on limall number of observation records Model output indicates surveys 
are most needed in and around Glaciei' National Pai'k, in the Seeley/Swan region, and 

portions of the South Fort: of the Flathead River drainage 


Oreohelix stn^osa 


G5/S5 


40 


39 


39 


Species IS clearly widespread 


Oreohelix stri^sa herry't 


G5T2/S1S2 


17 


8 


8 


Model based on limall number of observation records Model output indicates surveys 
are most needed m the Big and Little Snowy, Judith, Little Belt, Big Belt, eastern 
Elkhom and noilheni ed^e of the Beartooth Mountains 


Oreohelix stngosa depressa 


G5T5/SKR 


3 


3 


3 


Model based on exti^emely small numbei' of obseivati on records Species may best be 
detected dunng surveys for other species 


Oreohelix sid^nidis 


G5/S5 


49 


43 


43 


Species is clearly widespread Additional surveys in southwest Montana would seem 
likely to detect Iheir prei^ence 


Oreohelix ya\^apai 


G5/3NR 


4 








Not modeled due to limited numbei' of obseivation records 


OreoheUx yavapai mariae 


G4T1/S1 


5 


1 


"\ 


Model based on extremely lanall number of obi^ervation recordij Model output 
indicates suiveys are most needed in selected portions of the Bridger, Tobacco Root, 
Madison, and Gallatin Mountain Ranges 


Polygyrella polygyretia 


G3/S1S2 


13 


15 


15 


Model output indicates surveys are most needed in the Coeur d'Alene Mountains, the 
Bitteiroot Mountains of Minei'al County, portions of the Cabinet Mountains, and the 
area in and ai'ound Glaciei' National Parte 


Pi'istiloma ai'cticum 


G3G4ySMl 


1 


1 





Not modeled due to limited numbei' of obijervationrecordii 


Pimctism randolphi 


G4/SNR 


3 


3 


3 


Model based on extremely i^mall number of obi^ervation recordii Model output 
indicates surveys are most needed in the Coeur d'Alene Mountains and the Bitteiroot 
Mountains of Minei'al County 


P\q>illa hebes 


G5/3NR 


1 


1 





Not modeled due to limited numbei' of obseivation records 


Rndiodisci^ abietinn 


G4/S3S4 


79 


65 


57 


Model output indicates surveys are most needed in the Mission Mountains, the area in 
and around Glacier National Park, and the Bitteiroot di'ainages in Ravalli County 


Vallonia cyclophorella 


G5/SNR 


5 


5 


5 


Model based on small number of obseivation records Species may best be detected 
durinjj suiveys for other species 


Vallon ia grac the osta 


G5Q/3435 


27 


27 


27 


Species is clearly widespread Model output indicates suiveys are most needed in the 
eastern Elkhom Mountains, southern Big Belt Mountains, and near Bighom Canyon 


Vallonia pidchella 


G5/3NR 


1 


1 





Not modeled due to limited numbei' of obseivation records 


Veitigo crtstata 


G5/3NR 


1 


1 





Not modeled due to limited numbei' of obseivation records 


Veiti^o ^oiddii 


GS/SNR 


■h 


"\ 


1 


Very few records to base model on 


Veitigo modesta 


G5/34S5 


13 


13 


13 


Species IS clearly widespread Species likely went undetected a number of survey 
sites due to its small size and low detection probability 


Veiti^o ovata 


G5/3NR 


1 


1 





Not modeled due to limited numbei' of obseivation records 


Vitrina peilnada 


G5/S5 


111 


111 


111 


Species IS clearly widespread 


Zonitoides arboreiis 


G5/S5 


157 


151 


151 


Species is clearly widespread 


Zoogenetes harpa 


G5/3NR 


1 


1 





Not modeled due to limited numbei' of obseivation records 


SLUGS 1 



t 



Anon dtstinctus 


G5/3NA 


] 


1 





Not modeled Exotic species with limited observations 


Avionfasctatiis 


G5/3NA 


4 


4 


1 


Not modeled Exotio species with limited observations 


Anon mteiTnediiis 


G5/3NA 


4 


4 


■h 


Not modeled Exotic species with limited observations 


Anon nifiis 


G4G5/SNA 


5 


5 


5 


Not modeled Exotic species with limited observations 


Anon sisbfiiscus 


G5/3NA 


10 


10 


10 


Not modeled Exotic species with limited observations 


Deroceras laeve 


G5/SNA 


40 


40 


38 


Exotic species, but modeled because species appears to be widespread and its 

interactions with native species may become impoitant 


Deroceras pano rm itanum 


G5/3NA 


T 


"• 


■Y 


Not modeled Exotic species with limited observations 


Deroceras rettcidatitm 


G5/3NA 


7 


1 


7 


Not modeled Exotic species with limited observations 


HemphiiUa camehis 


G4/S1S2 


11 


11 


7 


Model based on small numbei' of obseivation records Model output indicates surveys 
are most needed in Lincoln and southwestern Flathead Counties 


HempkiUla danieisi 


G2G3/S1S2 


17 


12 


12 


Model output indicates surveys are most needed in the area in and around Glacier 
National Park and the Bitterroot Mountains in Minei'al and Ravalli Counties 


Kootenaia burkei 


G2/S1S2 


21 


21 


20 


Model output indicates surveys are most needed in the area in and around Glacier 
National Park and the Bitterroot Mountains in Mineral County, and the Seeley/Swan 
area 


Limax ma:^imiis 


G5/3NA 


11 


11 


10 


Not modeled Exotic species with limited observations 


Ma^ipelta mycopha^a 


G3/S2S3 


33 


30 


26 


Model output indicates surveys are most needed in the area in and around Glacier 
National Park, the Southfoiic of the Flathead drainage, the BitteiToot Mountains in 
Mineral County, and the Seeley/Swan ai'ea 


Prophysaon andersoni 


G5/S1S2 


7 


7 


6 


Model based on small numbei' of obseivation records Model output indicates surveys 
are most needed in the lower Clai'k Fork Rivei^ Bull River, and Lake Creek drainages 


Prophysaon humile 


G3/S2S3 


46 


45 


44 


Model output indicates suiveys ai'e needed in and around Glaciei' Nation Park, the 
Seeley/Swan area, the Mission Mountains, the east Cabinet Mountains, and the 
Bitterroot Mountains in Mineral County 


Udosarx lyrata 


G2/S1 


9 


6 


7 


Llodel output indicates suiveys ai'e needed in tlie Bitterroot Mountains in Minei'al and 
Ravalli Counties, and the area around the northern portions of the Flint Creek 
Mountains 


Zacoleus idahoensis 


G3G4/S2S3 


34 


27 


27 


Model output indicates suiveys ai'e needed in the Bitterroot Mountains in Minei'al 
County and on the western ed^e of the Mission Mountains 



Descriptions of Environmental Input Layers 

Input environmental layers consisted of 6 continuous (elevation, slope, curvature of land surface, an- 
nual precipitation, inaxiinuin July temperature, minimum Januaiy tempeiatiue) and 5 categorical (aspect, 
geology 1992 National Land Cover Data, soil temperature class, STATSGO soils) variables. These layers 
were tiled together, resampled, and converted to a statewide coverage of 90 meter grid cells in ASCII 
fomiatwitli 10,204 columns and 5,892 rows projected to North American Datum 1983 in Montana State 
Plane. Each source enviionmental layer is described below and, where appropriate, links to metadata are 
provided. 

Aspect (Categorical) 

Calculated using the Aspect fimction in AjxMap 9.2 Spatial Analyst from tlie 10-meter National Eleva- 
tion Dataset (NED) and resampled to 90 meters. See description of the elevation layer and the associated 
NED metadata link below. A brief summary of raster cell value descriptions for grid beaiing values fol- 
lows below. 

= Flat 

l=Norlh (337.5-22.5) 

2 = Noitheast (22.5-67.5) 

3 =East (67.5-112.5) 

4 = Southeast (112.5-157.5) 

5 = South (157.5-202.5) 

6 = Southwest (202.5-247.5) 

7 = West (247.5-292.5) 

8 = Northwest (292.5-337.5) 

Curvature of Land Surf ace (Continuous) 

Calculated using the Curvature fimction in ArcMap 9.2 Spatial Aialyst from tlie 90-meter resampled grid 
resulting from the 10-meter National Elevation Dataset (NED). See description of die elevation layer and 
the associated NED metadata link below. Values aie continuous from -50 (concave land surface) to +50 
(convex land surface). 

Elevation (Continuous) 

Tlie National Elevation Dataset (NED) is a 1/3 arc-second (lO-meter) raster grid of decimal meter values 
assembled by the U.S. Geological Sui"vey. Metadata on the 1/3 aic-second (10-nieter) NED is at: littp:// 
scamless.usgs.gov/prQducts/3aix.php 

Geology (Categorical) 

A polygonal coverage of suiUcial geology available in a mixture of 1:100.000 and 1:250,000 scales from 
tlie Montana State Geologic Mapping Program at die Montana Bureau of Mines and Geology. Metadata 
on the state geology map is at: http:/Av\vw.mbing. mtech.edu/gmr/gnir-statemap. asp 

Maximum July Temperature (Continuous) 

A polygonal coverage of estimated average maximum daily temperatures for July in degrees Fahrenheit, 
for tlie climatological period 1971-2000. Estimates are based on Parameter-elevation Regressions on 
Independent Slopes Model (PRISM) derived raster' data which uses known point temperature data and a 
digital elevation model (DEM) to generate gridded estimates of annual, monthly and event-based climatic 
paj'ameters. General infomiation on tlie underlying PRISM data and tlie source data itself can be down- 
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loaded from the Oiegon Climate Service website at: http:/Avww. ocs.orst.edu/prism/ . The Montana data 
reprojected to Montana State Plane and resampled to a resolution of 600 meters representing 33 tempera- 
ture ranges in degrees FaJii'enheit is available at: http://nris.int.gQv/nsdi/niis/tinax71_QQ.html 

Minimum January Temperature (Continuous) 

A polygonal coverage of estimated average minimum daily temperatures for January in degrees Fahren- 
heit, for tlie climatological period 1971-2000. Estimates are based on Parameter-elevation Regressions on 
hidependent Slopes Model (PRISM) derived raster data which uses known point temperature data and a 
digital elevation model (DEM) to generate gridded estimates of annual, monthly and event-based climatic 
parameters. General infomiation on tlie underlying PRISM data and tlie source data itself can be down- 
loaded fiom die Oiegon Climate Service website at: http ://www. ocs. or st . edu/prism/ . The Montana data 
reprojected to Montana State Plane and resampled to a resolution of 600 meters representing 33 tempera- 
ture ranges in degrees Falirenheit is available at: littp://nris.mt.gQv/nsdi/nris/tinin 7 l_OO.html 

National Landcover Data (Categorical) 

Tlie 1992 National Land Cover Data Set is based on 3Q-meter Landsat Tliematic Mapper imagery. Abrief 
summaiy of rastei' cell value descriptions follows below. Metadata on tliis layer can be found at: http:// 
nris.mt.gov/nsdi/nris/nlcdgrid.litml 

11 Open Water 

12 Perennial Ice/Snow 

21 Low Intensity Residential 

22 High Intensity Residential 

23 Commercial/Industrial/Transportation 

31 Bare Rock/Sand/Clay 

32 Quarries/Stiip Mines/Giavel Pits 

33 Transitional 

41 Deciduous Forest 

42 Evergreen Forest 

43 Mixed Forest 
51 Shi'ubland 

61 Oichards/Vineyards/Otlier 

71 Grasslands/Herbaceous 

81 Pasture/Hay 

82 Row Q'ops 

83 Small Grains 

84 Fallow 

85 Urban/Recreational Grasses 

91 Woody Wetlands 

92 Emergent Herbaceous Wetlands 

Annual Precipitation (Continuous) 

A polygonal coverage of annual precipitation in inches for the climatological period 1961-1990 based 
on source data fiom the National Weatlier Seivice Cooperative stations. Natural Resources Conservation 
Service SNOTEL stations, and local networks. Metadata on tliis layer can be found at: http ://ni is.nit.gov/ 
nsdi/nris/precip .htm I 
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STATSGO Soils (Categorical) 

State Soil Geographic data (STATSGO) is a polygonal coverage of general soil associations developed by 
the National Cooperative Soil Sui'vey Tlie soil maps for STATSGO aie compiled by generalizing more 
detailed soil survey maps. Map unit composition for a STATSGO map is determined by transecting or 
sampling areas and expanding tlie data statistically to chaiacterize tlie whole map unit. Tlierefore, map 
units depict the dominant soils making up the landscape and often contain dissiinilai' soil types. The ap- 
proximate minimum aiea delineated is 625 hectai'es (1,544 acres). Background infonnation and metadata 
on STATSGO is available at: http://nris.mt.gov/nsdi/statsgo.pdf . http://nris.mt.gov/nsdi/nris/SS19.html , 
and http://dbw^v^v.essc.psu.edu/doc/5tatsgo/statsao info.html . Definitions for tlie 694 map units used in 
tlie input can be downloaded as .dbf files along witli a STATSGO shapefile for Montana at: littp://nris. 
iiit.sQv/gis/sisdatalib/sisDataList.aspx?datasroup=state\^ude-resionalfeeaichTenns=statsgo 

Slope (Continuous) 

Percent slope (rise over run multiplied by 100) calculated using the Slope function in Ai'cMap 9.2 Spatial 
Analyst fiom tlie 10-meter National Elevation Dataset (NED) and resampled to 90 meters. See descrip- 
tion of the elevation layer and tlie associated NED metadata link above. 

Soil Temperature Regime (Categorical) 

A generalized polygonal coverage of soil temperature regimes. Abrief summary of raster cell value 
descriptions follows below. A glossary of relevant soil tenninology can be found at: https://www.soils. 
org/sssasloss/ . Metadata on tliis layer can be found at: littp://soils.usda Qov/use/worldsoils/mapindex/str. 
html 

1 Qyic/Udic 

2 Frigid/Udic 

3 Frigid/Typic Ustic 

4 Ciyic/Typic Ustic 

5 Frigid/Aridic Ustic 

6 Frigid Aquic 

7 Frigid/Typic Xeric 

8 Watei' 

9 Ciyic/Typic Xeric 

10 Ciyic/Ai'idic Ustic 

11 Mesic/Ustic Aridic 

12 Ciyic/Udic Ustic 
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Example Interpretation of Model Output 

Tlie main map output tliat we show is in thefonn of a hot-to-cold color map indicating tlie suitability of 
each grid cell as a function of the environmental vaiiables at tliat grid cell (Figure El). A high (hot) value 
of the function at a paiticulai' grid cell indicates that the grid cell is predicted to have suitable conditions 
for tliat species. Thus, grid cells having values of 100 represent the most suitable habitat (hot or best fit 
to the positive data) while tliose witli values closer to represent less suitable habitat (cold or worse fit 
to the positive data). Note tliat this output is not inteipretable as the probability of a species' occupancy. 
Rather it is a continuous cumulative function with tlie cumulative value of a cell indicating tlie percent- 
age of grid cells tliat have raw outputs fiom tlie Maxent exponential model tliat are the value of tliat cell 
or smaller Perhaps tlie best interpretation of the cumulative output is in tenns of tlie expected omission 
eiTor rate (tlie rate at which obseiTations occur in areas where tlie model does not predict tliem). Thus 
a cumulative value threshold of X would be expected to have an X% omission eiror rate. Because most 
ten'eshial mollusk species have veiy little data and because tlie main goals of tliis modeling effort were 
to identify high probability ai'eas tliat cuiTently lack observation infomiation for individual species, we 
truncated tlie continuous fiinction at a tlireshold that balances tlie probability of omission and commission 
error rates and show only two levels of model output (orange and red) tliat aie predicted to have tlie most 
suitable habitat conditions for the species (Figure El). 




Figure EL Predicted distribution model for Smoky Taddropper {Prophysaon hwnile). Hotter colors indicate areas that are 
predicted to have more suitable habitat for the species. Black dots are positive data i^edto build the model. Gray dots are loca- 
tions where a timed visual encounter siu-vey for terrestrial mollusks was performed during the 2005-2007 field seasons. Land- 
ownership, a shaded relief map, and county lines are included for reference. 
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Figure E2. An evaluation of omission error rates for training data as afimction of the <:itmidatwe threshold and overall predicted 
area. Note that all data available for terrestrial mollasks was used to tram the models so there was insufficient data to indepen- 
dently test the models. Ihe red line indicates the overall fraction of the map area fitting each value of the cwnidative threshold. 
Thus, for the Smofy Taildroppei; about 70% of Montana has a cumulative threshold near zero (i.e. no suitable habitat). Ihe light 
blue line is the predicted omission rate for each cwnidative threshold. Notice that the actual omission rate of training samples 
(dark blue line) is lower than the predicted omission rate with no omissions occurring until a cianidative threshold value of 3 and 
an overall fractional map area of 0.025. 
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Sensitivity vs. Specificity for Smoky Taildropper 
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Figure E3. A Receiver Operating Characteristic (ROC) cwve evaluating the overall predictive power of the model with the Area 
Under the Curve (AUCj. Ihe AUC value for the Smoky Taildropper model of 0.996 indicates that when iivo random locations 
are chosen the model has a 99.6% prohahilit}' of assigning a higher cwnidative threshold value to the location with more suit- 
able habitat. The light blue line indicates how a neidral or random model would perform (i.e., it only has a 50% probability of 
assigning a higher cuimdative threshold value to a random location with more suitable habitat than a random location with less 
suitable habitat). Ihe fitither toward the top left of the graph the red training data line is. the better the model is at predicting the 
presences contained in the training data. Sensitivity (plotted on they-axis) is the prop oition of positive locations that were cor- 
rectly classified by the model. Sensitivity is also known as the tnie positive rate and can be thought of as the degree of absence 
of omission errors. Specificity is the proportion of random locations chosen from the background (these pseudo-absences are 
used instead of tnie negative locations) that were correctly classified by the model as negative. One minus the Specificily (plotted 
on thex-axis) is known as the false positive rate and represents the commission error rate. Ihi^, the Smoky Taildropper model 
below has a high absence of omission errors and a low occurrence of commission errors. 
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Figure E4. Response cuivesfor individual environmental variables showing how the logistic prediction changes as each environ- 
mental variable is varied while all other environ/jjental variables are held constant at their average sample values. The value 
on they-axis is predicted probability^ of suitable conditions as given by the logistic fonmda P(x) = exp(cl ^fl(x) + c2 *fi{x} + 
c3 *f3(x) )/ Z. Note that if any of the environmental variables are correlated, the maigmal response cwves can be misleading 
(e.g., two highly correlated variables with opposite response cwves could effectively cancel each other out). Value definitions 
for categorical variables are either shown in the descriptions of the environmental variables or can be accessed via the metadata 
links provided. 
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Figure E4. (Contimied). Response ciwves for individual environmental variables showing how the logistic prediction changes as 
each environmental variable is varied while all other environmental variables are held constant at their average sample values. 
Ihe value on they-axis is predicted probability of suitable conditions as given by the logistic fonnida P(x) = exp(cl '^'fl{y.) + c2 
*f2{xj + c3 ^y^fx) )/Z. Note that if any of the environmental variables are correlated, the marginal response ciuves can be 
misleading (e.g., two highly correlated variables with opposite response curves could effectively cancel each other oiU:}- Value 
definitions for categorical variables are either shown in the descriptions of the environmental variables or can be accessed via 
the metadata links provided. 
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Figure E5. Jackknife chart showing the relative impoitance of environmental variables as afimction of the change m "gain" 
(the log of the number of grid cells minus the average of the negative log probabilities of the sample locations) resulting from the 

exclusion or sole inchision of the environmental variable in the model. Variables with the highest training gain resulting from 
sole inclusion of those variables (dark blue bars) are the best individual variables at describing suitable habitat for the species. 
Variables with the greatest reduction in training gain resulting from their exclusion {light blue bars) contain information on the 
species habitat u^e that is not present in other variables. The red bar indicates the maximum gain achieved with inclusion of all 
variables. 
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